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Evolution of the scale factor a{t) in Friedmann models (those with zero pressure and a constant 
cosmological term A) is well understood, and elegantly summarized in the review of Felten and 
Isaacman [Rev. Mod. Phys. 58, 689 (1986)]. Developments in particle physics and inflationary 
theory, however, increasingly indicate that A ought to be treated as a dynamical quantity. We 
revisit the evolution of the scale factor with a variable A-term, and also generalize the treatment to 
QQ ' include nonzero pressure. New solutions are obtained and evaluated using a variety of observational 

I criteria. Existing arguments for the inevitability of a big bang (ie., an initial state with a = Q) are 

^\ . substantially weakened, and can be evaded in some cases with Aq (the present value of A) well below 

current experimental limits. 
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I. INTRODUCTION 



The behavior of the cosmological scale factor a{t) in solutions of Einstein's field equations with the Robertson- 
Walker line element has been the subject of numerous studies. Textbook presentations tend to focus on models in 
which pressure p is zero and there is no cosmological term (A = 0). Some treatments include the Friedmann models, 
in which p = but A 7^ Much less attention has been directed at the more general Lemaitre models, in which 

• pressure p is given in terms of density p by an equation of state p = p{p), and A 7^ 

■ The possibility of a nonzero A-term, in particular, has resurfaced lately in connection with the age problem [ p^ . If 
00 A is large enough, in fact, the age of the Universe (defined in the standard model as the time elapsed since a = 0) 

■ can in principle become infinite, as a{t) never drops below a nonzero minimum value a* in the past direction. The 
(— I , existence of such "big bangless" oscillating models has been recognized for over sixty years ||l2| . They have however 

i' been dismissed as unphysical on the grounds that the required values of A are incompatible with observation [p^[-p^ . 
' We will return to the observational constraints later on; merely noting here that the above arguments, along 
^ • with nearly all existing astrophysical work on the cosmological term, operate on the assumption that A = constant. 
' Quantum field theorists and others, by contrast, have been treating the cosmological term as a dynamical quantity for 
; thirty years [^^-^. Anything which contributes to the energy density p„ of the vacuum behaves like a cosmological 
^ term via A„ — SnGpy. Many potential sources of fluctuating vacuum energy have now been identified, including 
• I— I , scalar fields , tensor fields |^,^ , nonlocal effects [|7| , wormholes |28) , inflationary mechanisms ^^^d even 

' cosmological perturbations . Each of these has been claimed to give rise to a negative energy density which grows 
^ , with time, tending to cancel out any pre-existing positive cosmological term and drive the net value of A toward 
■ - - I zero. Processes of this kind are among the most promising ways to resolve the longstanding cosmological "constant" 
problem (see for review). 

The purpose of the present paper is to re-examine the evolution of the scale factor when A 7^ constant. This has not 
yet been done in a systematic way. We also expand on most earlier treatments by considering a fairly general equation 
of state for ordinary matter rather than restricting ourselves to the pressure-free Friedmann models. These two 
generalizations lead to qualitatively as well as quantitatively new behavior for a{t), and hence for related phenomena 
such as the age of the Universe, the appearance of gravitational lenses, and the redshifts of distant astronomical 
objects. They can also allow one to circumvent the abovementioned arguments against oscillating models. 

The remainder of the paper is organized as follows: the required assumptions, definitions and dynamical equations 
are assembled in §||, and applied in §§[!! - Vl to models in which A varies as a function of time t, the scale factor 



a, the Hubble parameter H and the deceleration parameter q respectively. In each section we obtain analytic or 
numerical solutions for the scale factor, paying particular attention to the question of the initial singularity, and 



discuss observational constraints where appropriate. Conclusions are summarized in §VII 



II. VARIABLE-A COSMOLOGY 
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A. The Cosmological "Constant" 



To begin, we recall why the cosmological term has often been treated as a constant of nature. The Einstein field 
equations read: 

G^. + A.g^, = 8^GT,,,, (2.1) 
where G^i^ = — is the Einstein tensor and T^^ is the energy-momentum tensor of matter (we assume units 



such that c = 1). Taking the covariant divergence of eq. (2.1), recalling that the vanishing covariant divergence of 



the Einstein tensor is guaranteed by the Bianchi identities, and assuming that the energy-momentum tensor satisfies 
the conservation law V^T^jy = 0, it follows that the covariant divergence of A^^^ must vanish also, and hence that 
A = constant. This argument, which situates A firmly on the left-hand side of the field equations, constitutes a 
"geometrical interpretation" of the cosmological term. 
More rec entl y, it has become increasingly common (see, eg., p2[) to move the cosmological term to the right-hand 



side of eq. (|2.lD 



G^v — SnGT^i,, (2.2) 

where: 

ottG 

that is, to interpret A as part of the matter content of the universe, rather than a purely geometrical entity. Once 
this is done, there are no a priori reasons why A should not vary — as long as the effective energy-momentum tensor 
T^y satisfies energy conservation Q. 



B. Dynamical Equations 



We will make the usual assumptions: a homogeneous and isotropic Universe (ie., Rob ertso n- Walker line element) 
and perfect-fluid- like ordinary matter with pressure p and energy density p. Eq. (2.3) then implies that the 
effective energy- momentum tensor also has the perfect fluid form, with effective pressure p = p — A/SttG and energy 
density p = p + A/SttG The field equations (2.1) and law of energy conservation then read: 



SttG 



A 2 

—pa- + 3 a 



d 

da 



3 p 



k 

A 
8^ 



For the equation of state we take: 



P = (7 - 



(2.4) 
(2.5) 

(2.6) 



with 7 = constant. Previous analytic and numerical studies of the evolution of the scale factor have tended to focus 



on the dust-like case 7 = 1 and occasionally also the radiation- like case 7 = 4/3 j6|-|^ 

ourselves to these values, as many other possibilities have been considered in the literature p7 
Substituting eq. (p^) into eq. (|2.5|), we find 



a^f dA 



We will not restrict 



(2.7) 



When A = constant, eq. (2.7) reverts to the well-known result that density p scales as a^"^ in a pressure-free universe 
(7 = 1) and a~* in a radi ation-dominated universe (7 = 4/3). 

Differentiating eq. (2^) with respect to time and inserting eq. (p/7|), we obtain: 



SttG 



1-1 



pa 



A 

3"- 



(2.8) 



Eq. (2.8) shows that a positive density p acts to decelerate the expansion, as expected — but only when 7 > 2/3. If, 



on the other hand, the cosmological fluid is such that 7 < 2/3, then its density can actually accelerate the expansion. 
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a > (this is the phenomenon commonly known as inflation). Eq. (2.8) also confirms that a positive cosmological 
term contributes positively to acceleration, "propping up" the scale factor against the deceleration caused by the 
matter term p. A negative A-term, on the other hand, acts in the opposite direction and brings about recollapse more 

quickly. 

Combining eqs. (2.4) and (2.8), we can eliminate p to obtain a differential equation for the scale factor in terms of 
the cosmological term alone: 



(2.9) 



This differential equation governs the behavior of the scale factor in the presence of a cosmological term A, whether 
or not the latter is constant. 



C. Phenomenological A-Decay Laws 



The abovementioned sources of negative vacuum energy |23|-|3C|] do not, in general, lend themselves to simple 
expressions for A in terms of t, a, H or q. There are some exceptions, including scalar field-based ||l^,|2^,^,^ and 
other theories ||4^ , |53| , |7^ -f7^ in which analytic decay laws are derived from modified versions of the Einstein action. In 
most such papers, however, no exact solution for A is obtained; the intent is merely to demonstrate that decay (and 
preferably near-cancellation) of the effective cosmological term is possible in principle. 

In a complementary approach, a number of authors have constructed models of a more phenomenological character, 
in which specific decay laws are postulated for A within general relativity. These theories are incomplete to the 
extent that they do not include explicit physical mechanisms to govern the exchange of energy between the shrinking 
cosmological term and other forms of matter |Q. In some models this issue is not addressed at all; in others the 
energy is assumed to be channelled into production of baryonic matter and/or radiation. The former case can be 
constrained by observations of the diffuse gamma-ray background, since (assuming the decay process does not violate 



baryon number) one might expect equal amounts of matter and antimatter to be formed 1 47 59[| . The latter can be 
constrained by nucleosynthesis arguments cosmic microwave background (CMB) anisotropics [ ^Tp^JS^ , the 

absolute CMB intensity [^0|j59| ] , and thermodynamical considerations |5^j6^,^ . 

We take the point of view here that simple A-decay scenarios are worth examining, irrespective of whether they come 
from extended theories of gravity or phenomenological considerations, for at least four reasons: (1) they have been 
shown to address a number of pressing problems in cosmology (p8[| - |]77[|); (2) many are independently motivated, 
eg., by dimensional arguments, or as limiting cases of more complicated theories; (3) most are simple enough that 
meaningful conclusions can be drawn about their viability; and (4) successful implementation would point toward the 
eventual Lagrangian formulation of a more complete theory. For convenience, we have collected together the most 
common decay laws from the literature and listed them in Table || (by chronological order of appearance) . 

In the remainder of the paper, we focus on power-law functions of one parameter. Our discussion is however more 
general than most of those noted in Table |, because we do not fix values of the exponents a priori. In particular, we 
consider decay laws of the following four kinds: 

A = At-^ (2.10) 
A = i3a-™ (2.11) 
A = Ci/" (2.12) 
K^Vq\ (2.13) 

where y^, S, C, I?, ^, to, rt and r serve as adjustable constants. These four variable-A scenarios and their cosmological 
consequences are explored in §111, §|v|, §M and §VI respectively. 



D. Definitions 

We conclude §^ by introducing the terms and definitions which will be needed to connect our solutions to obser- 
vation. Chief among these are the energy densities of ordinary matter and the cosmological term, expressed in units 
of the critical density: 



A ^ iHl^ 

Peril SHq SttG 



' ^ = ^7772 ' Pcru = ^. (2.14) 
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We will be particularly interested in the values of these parameters at the present time (subscript "0" ) : 

These will constitute our primary free parameters throughout the following sections. 

The usef ulness of the quantity Aq (sometim es d enoted VIa in the literature) is highlighted by using the defini- 



tions (2.14) to rewrite the Lemaitre equation ( |2.4| ) in the form = Hq{Q + A) — k/a^. At the present time this 
implies: 

k = alH^{no + Xo-l). (2.16) 



From eq. ( ^.16 ) it is clear that: 



f7o + Ao>l =^ A:>0 
f7o + Ao = 1 =^ fc = 

l}o + Ao<l =^ fc<0. (2.17) 

Most cosmologists implicitly choose units for uq such that the value of fc is normalized to either or ±1. We will 
follow Felten and Isaacman in refraining from this, because it is more convenient for our purposes to choose units 
such that flo = 1. 

In place of and A, some authors prefer to use the quantities a and q, defined by: 



AttGo da /„ „x 

, (2.18) 



The present values of these parameters are fixed with the help of eqs. (2.8) and(2.14) 



ao = r!o/2 , go = (37/2 - l)r!o - Ao. (2.19) 



Negative values of qo can be obtained for large (positive) Aq, or small 7. Eq. ( ^.19 ) yields the standard expression j|] 
for qo when 7=1. 

III. A AS A FUNCTION OF TIME 
A. Interpretation of the Time Co-ordinate 

We begin with the oldest, and probably the most straightforward implementation of the variable- A idea; namely. 



one in which the cosmological term is a simple power-law function of time, as set out in eq. (2.10) 



A^At-^, (3.1) 

where, for later convenience, we measure time in units of Hubble times (r = HqI). The case £ — 2 has previously 
been considered by several authors. Each, however, has imposed supplementary conditions. Endo, Fukui and others 
p9| , p4t operate in the context of a modified Brans-Dicke theory. Canuto et al assume invariance under changes of 
scale, while Lau adopts the Dirac large-number hypothesis (with a time-varying gravitational "constant" G) from 
the outset. Berman requires that the density p of ordinary matter also scale as r~^, and that the deceleration 
parameter g be a constant. Beesham ^o|| restricts his treatment to Bianchi Type I models with variable G. Lopez 
and Nanopoulos ]4l| ] take A to have the same dependence on the scale factor (A cx a~^) as on time, for late times 
at least. (These latter authors also make the important claim that a A-decay ansatz of this kind could follow from 
certain versions of string theory.) The question of the initial si ngul arity is not addressed in any of these papers. In 



this paper, we examine the properties of models with the form (3.1) in a more comprehensive way. 

Several conceptual issues should be dealt with before we proceed. Firstly, since we are interested in oscillating as 



well as traditional "big bang" models, we are obliged to broaden the conventional definition of cosmic time |78|, in 
which it is set to zero at the moment when a = 0. In those cases where the scale factor never vanishes, we choose 
here to measure r instead from the most recent moment when da/dr = (the "big bounce"). In either case we refer 
to this as the "initial moment" and denote it by t = t, = 0. (Here and elsewhere, the subscript "*" will indicate 
quantities taken at this time.) 
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Secondly, eq. (3.1) implies that A — !■ cx) in the limit r — > 0, which may not be realistic. A may go to infinity at some 
other time (tqo, say) than the initial moment. Or, more likely, its divergent behavior may be cut off at some critical 
temperature (at time Tc, say) by a pha se tr ansition or similar mechanism, above which A is effectively constant [ plf . 
So a more plausible formulation of eq. (3.1) might take the form: 



A 



Ae 



when T < Tc 
when T > Tr 



(3.2) 



where continuity across the cutoff time Tc implies that Ac = A{tc — Too)~^- However, the fact that the dynamical 
equations of contain no explic it ti me-dependence means that we c an shift our time co-ordinate via (r — Tqo) — *■ t 
with impunity. The decay law ( |3.2| ) then reverts to the form (3.1), for all r > Tc at least. In practice we will 
assume that eq. (3.1) holds for all times of in terest. For earlier times the standard (Lcmaitre) solution applies, with 
A = Ac = constant. If Tc~ Too, then eq. (|3.lD holds all the way to r w 0. 



B. Riccati's Equation 

We are now in a position to study the evolution of the scale factor. Switching dependent variables from a to H: 

H^^-=hJ^), (3.3) 
a \aaT I 



we obtain for the differential equation (2.9): 

dH _ f-i-f 

d7 ^ ynfo 



^ ) A 



V2i/o 



1 - 



37 

2 ) Hqu^' 



(3.4) 



If we restrict ourselves to spatially flat universes (fc = 0), then the last term drops off, leaving a special case of Riccati's 
equation: 



fl J-f 

— =7'(r) H^ + QiT)H + n{T), 
ar 



(3.5) 



where V = —3j/2Ho, Q = and TZ(t) = {-^ /2Hq)K{t). We will adopt this restriction for the remainder of §111. This 
has also been imposed in most of the special cases studied so far [ pi)p^j3[^ , [Il[ | . 

Solving Riccati's equation by standard methods |f79|, we change dependent variables from H to x via: 



^_ I dx _ f2Ha\ dx 
Vx dr \ 3j J X dr^ 



(3.6) 



whereupon eq. (3.4) takes the form [putting k = and inserting eq. (3.1) for A(t)] 



with: 



T — - ax = 0, 
dr'^ 



4i/2- 



(3.7) 



(3.8) 



This is now linear, as desired. We will solve for x{t) in the cases ^ = 1,2, 3, 4. 

Once x{t) is found, the Hubble parameter H{t) follows imme diat ely fro m th e definition ( |3.6| ). Moreover, the scale 
factor air) is also known, as may be verified by comparing eqs. (3.3) and (|3.6|) to yield: 



a(r) = [x(r)]2/37. 



(3.9) 



The con stant a given by eq. (3^) can be fixed in terms of observational quantities as follows. Applying the decay 
law (3.1) to the pres ent epoch t = rg, and using the definition (2.15) of Aq, we find that A = ZHqXqTq, which can be 
substituted into eq. (3_^) to yield: 



a = (37/2)^AoTo^. 



(3.10) 



With the restriction k = 0, eq. (2.17) implies that VIq is given by (1 — Aq) throughout §111 
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C. The Case £ = 1 



We now proceed to the first case. The differential equation (3/7) reads: 



dT^ 



ax = 0, 



where a is given by eq. ( 3.10 ): 

a = (37/2)2Aoro. 

Following standard techniques [fzof , we switch independent variables from t to z = 2y/—aT, whereupon: 



d X dx 2 

-r-^ — z - — \- z x — 0. 

dz'^ dz 



(3.11) 



(3.12) 



(3.13) 



This is transformable to Bessel's equation, with general solution x{z) = cizJi(z) + C2zYi(z), where Ji{z) and Yi{z) 
are Bessel and Neumann functions of order one. Eq. (|3.9|) then gives for the scale factor: 



a(r) 



- ^1/37 



[ciJi(z)+C2Yi(z)] 



(3.14) 



where we have absorbed a factor of 2y^—a into ci, C2. The Hubble parameter is found by putting x{z) into eq. (t3.6|): 

ciJoiz) + C2Yo{z) 



H{t) = Ho 



To 



ci Ji(z) + 02^1(2;) 



(3.15) 



where Jo{z) and Yo{z) are Bessel and Neumann functions of order zero. We note from the definition ( 3.12 ) that z{t), 
and hence a(T) and H{t), can only be real (for positive times) if Ao < 0, which would imply a negative cosmological 
constant. While this possibility has been considered in some contexts [|7^,Q, we will see shortly that it leads in the 
present theory to unrealistically short ages for the Universe. Therefore the case £ = 1 is probably not realized in 
nature. 

We can impose the following boundary conditions at the present epoch: 



a(To) = ao = 1 , H{tq) = Hq. 
Substituting these expressions into eqs. ( |3ll ) and dslsl ), it is straightforward to solve for ci and C2: 

^^Yoizo)-Yiizo) 



ci = 



C2 



V-AqTo [Ji{zo)Yo{zo) - Jq{zq)Yi{zo)] 

-V-Aq Jo(zo) + Ji{zo) 

V-AqTo [Ji{zo)Yo{za) - 70(2:0)^1(20)] ' 



where: 



20 = S"fTQy/-Xo. 



(3.16) 

(3.17) 
(3.18) 

(3.19) 



To keep a^r) finite at t = 0, we require C2 = 0, since Yi{z) diverges logarithmically at 2 = 0. This constitutes our 
third boundary condition. In conjunction with eq. ( 3.1§| ) it implies that Ji(zo) — V— Ao ^0( 20) = 0. This equation 
may be solved numerically fo r rn as a function of Ao, with the help of the definition (3.19). The results can then 



be substituted back into eq. ( 3.17 ) to fix the value of ci. With ci and C2 both known, a(r) and H(t) are given by 
eqs. (3.14) and (3.15) respectively. 

The evolution of the scale factor for this case is illustrated in Fig. |l| for various values of Ao, assuming 7 = 1 (after 
Felten and Isaacman Q). In particular, we have followed these authors in plotting the scale factor a as a function of 
(r — To), rather than r for each curve. This has the effect of shifting all the curves so that they intersect at (0, 1), 
which marks the present epoch. [Recall that we have chosen units such that ao = 1, eq. (2.16).] 

We have plotted for four Hubble times into the future, and one Hubble time into the past. It may be seen that 
the Ao = curve (solid line) intersects the time axis at (r — tq) = —2/3, which confirms the well-known rule that the 
age of a flat universe with no cosmological constant is tq — 2/3. The models with Ao < (dashed lines) are younger 
than this, which considerably diminishes their attractiveness. The Ao = — 1 model, for instance, has To = 0.48, while 
the Ao = —3 model has tq = 0.35. If we use a current widely-accepted value of Hq = 73± 10 km s^^ Mpc^^ pi| for 
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the Hubble parameter, then (recaUing that t = t/Hq) we see that the age of the Universe in these models can be no 
more than 7.4 and 5.4 billion years old respectively. This conflicts badly with estimated globular cluster ages, which 
are thought to be at least 9.6 billion years old in some cases |82[. The situation improves slightly if one switches to 
the lower values for Ho which are reported by some authors ||83|] 7^ If Hq = 55 ± 10 km s~^ Mpc^^, then the maximum 
possible age for these two models increases to 10.4 and 7.6 billion years respectively. One can safely rule out models 
with Aq < — 1 on this basis, while models in the range — 1 < Ao < remain marginally viable at best. 

Fig. |l| tells us that, while flat models with Aq = will continue to expand indefinitely as us ual, those with negative 
values of Aq will experience eventual recollapse. This can be understood by looking at eq. (2.8), which shows that 
there are two contributions to the deceleration: one (which goes as —pa) due to ordinary matter and the other (which 
goes as +Aa) due to the negative cosmological constant. Because the density p of matter rapidly thins out with 
expansion, the first contribution alone is not enough to close the Universe when fc = 0. The second contribution, 
however, is diluted much more slowly (A drops off as only in this case), and is therefore sufficient to turn the 
expansion around eventually, no matter how small its value at the present time. Thus, models with Aq = — 1 and 
Ao = —3 encounter the "big crunch" after only 2.94 and 1.21 Hubble times respectively. 



D. The Case £ = 2 



For this case, eq. ( p. 71) takes the form: 



with a given by eq. ( |3.10|) as follows: 



2 d'^x 

T -—r - ax = [), 



= (37/2)^ Aor2. 



(3.20) 



(3.21) 



This is a special case of Eulcr's differential equation. Applying standard methods [79| , we switch independent variables 
via 2/ = In r to recast eq. ( 3.2C ) in the form: 



d^x dx 
dy^ dy 



— a X ^ 0. 



(3.22) 



This now has constant coefficients, as desired. There are three subcases, according as Aq is greater than, equal to, or 
less than — 1/(37To)^. 



1. The Subcase Ao > -l/(37ro)^ 



Since we exp ect on observational grounds that Aq is probably positive, this is the most physical of the subcases. 
Solution of eq. (3.22) for x{y) and hence x{t) is straightforward. The scale factor and Hubble parameter are found 



from eqs. (3^) and ( |3.6| ) to be: 



i(r) = t1/37 (cir'"" + 



C2T 



,2/37 



2Ho 

37 



rriiCiT 



m2C2T" 



■ (ciT™" + C2T-™0) 



(3.23) 
(3.24) 



where toq = |-\/l + (37To)^Ao. It is clear that a(T) diverges at t ^ fo r Aq > unless C2 = 0, which we consequently 
adopt as a boundary condition (as in the £ = 1 case). Eqs. ( ^.23 ) and (3.24) simplify to: 



a(T) = (cir"! )2/3'^ 



H{r) 



2Ho 

37 



(3.25) 
(3.26) 



where mi = 1/2 + nio- We then apply the boundary conditions ( 3.16 ) aX t — tq. Eq. ( 3.26| ) with H{tq) — Hq implies 
that: 



To = 2mi/37. 



(3.27) 
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Sub stitut ing this result into eq. ( 3.25| ), with a(To) = 1, we find that ci = {3j/2mi)"^'^ , which can be put back into 
eq. (3.2£) to yield the following expression for a(T): 



2m 1/37 



(3.28) 



The scale factor expands as a simple power-law function of time. This is consistent with previous special cases 
obtained for £ = 2: Endo and Fukui's a(T) cx r^^/Hr^-i} Berman and Som's a(T) cx r^/™ Q, Lau and Beesham's 
Herman's a(r) cx t^^^ and Lopez and Nanopoulos' a(T) c>c r [^ . 
In conjimction with the definitions of mg and mi, eq. ( [3.27 ) fixes the age of the Universe at: 



To 



37(1 -Ao)' 



(3.29) 



from which we draw two conclusions: firstly, that all models satisfying the boundary conditions obey Aq < 1; and 
secondly, that the age of the Universe tq ^ 00 as Aq ^ 1. The initial singularity can thus be pushed back arbitrarily far 
into the past. We also find a lower limit on the age of the Universe in these models by noting that mi = (toq + 1/2) > 
1/2. Inserting this into eq. (3.27) produces the result: 



To > 1/37. 



(3.30) 



In other words, assuming zero pressure (7 = 1), all models in this case have survived for at least one-third of a Hubble 
time. Adopting a recent observational upper limit of 83 km s~^ Mpc ~^ ||8l|] , this implies a minimum age of at least 
3.9 bilhon years. Finally, putting eq. ( [3.30 ) back into the expression (3.2£) for tq, we find that Aq > —1, which fixes 
the critical value of Aq separating this subcase from the other two. 



2. The Subcase Aq = -l/(37ro)^ 



Solution of eq. ( 3.22 ) is also straightforward, and one finds from eqs. (3_^) and (3^) the following general expressions 
for a(T) and H{t): 



i(r) =T^/37 (c3-}-c4lnr)2/3^ 



H{t) = 



2Ho 
37 



(c3 + 2C4) + C4 Inr 



2t(c3 -I- C4 Inr) 



where C3, C4 are arbitrary constants. To keep a(T) finite at r = we require C4 = 0, so that: 

air) = (c3V^)2/3^ 



Inserting H{tq) — Hq into eq. (3.34), we find for the age of the Universe in this model: 

To = 1/37, 



(3.31) 
(3.32) 

(3.33) 
(3.34) 

(3.35) 



which is exactly the limit Ao ^ — 1 in eq. ( 3.29| ). Eq. ( |3.35|) also corresponds to the lower limit allowed by eq. ( ^.30 ) 
as one might expect. 

Substituting the age ( 3.35 ) into eq. ( 3.33 ), meanwhile, and imposing a(To) = 1 as usual, we find that C3 — 



Inserted back into eq. (3.33), this yields 



air) 



1/37 



(3.36) 



which joins smoothly onto the result (3.28) from the previous subcase 
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3. The Subcase Ao < -l/(37ro)^ 



Solution of eq. ( 3.22 ) is also straightforward and leads to the following general expressions for a^r) and H(t): 



air) 
H{t) 



rl/37 



[c5 sin(m3 In r) + ce cos(to3 In r)] 



2/37 



2Hq ( (c5 — 2m3C6) sin(m3 Inr) + (cg + 2TO3C5) cos(m3 Inr) 



37 



2t[c5 sin(m3 Inr) + cg cos(m3 Inr)] 



(3.37) 
(3.38) 



where = ^ (37ro)^Ao — 1 and C5,C6 are arbitrary constants. Application of the boundary conditions ( 3.16| ) 
gives C5 and cg in terms of tq : 



C5 = —= sin(m3 In To) H 



1 



eg = ^= cos(m3 Inro) 

Vn) W3 



37 



1 



37 



cos(m3 Inro) 
sin(m3 Inro). 



(3.39) 
(3.40) 



As usual, we need a third boundary condition to fix the value of ro. Unlike the previous two subcases, we cannot 
keep a(r) finite at r = by setting one of C5,cg to zero. Instead we have adopted a numerical approach, searching 
iteratively for the value of ro consistent with the boundary conditions (ie., with the requ ireme nt th at either a(r) or 
H{t) go smoothly to zero as r — > 0). The values of C5 and cg then follow from cqs. (3.4C) and (3.4C) respectively. 

The evolution of the scale factor for this case is illustrated in Fig. ^, which has exactly the same format as Fig. 
except that we have plotted for three Hubble times into the past instead of one. Fig. exhibits a richer variety of 
solutions than Fig. |l]. The most noticeable difference is the existence of solutions for positive Ao (short-dashed lines). 
Of particular interest is the limiting case Ao = 1, which only approaches a = asymptotically as r ^ — 00. This case 
is not very realistic, however, as it has zero density (since fio = 1 — Ao). It is in fact the empty de Sitter model. The 
same solution is found for these values of Ao and CIq in conventional Lemaitre cosmology with A = constant Q . 

Fig. |2| therefore shows that we cannot avoid the big bang in a theory with A oc r~^ and fe = 0. We can, however, 
significantly extend the age of the Universe. Suppose we choose Ao = 0.5, for example; a value compatible with the 
tightest observational bounds thus far |84|. Fig. ^ shows that this model would have come into being 1.33 Hubble 
times ago [see also eq. (3.29) above]. Even if we adopt the upper limit of 83 km s~^ Mpc^^ for Hq [81|]. this translates 
into an age of 15.7 billion years — more than enough time for the oldest globular clusters to form |p3]. By way of 
comparison, a constant-A model with Ao = f^o — 0.5 has an age of only 0.83 Hubble times. 

Flat models with no cosmological constant are represented in Fig. ||by the curve labelled Ao = (solid line). As in 
the £ = 1 case, these have an age of 2/3 Hubble times. Models with a negative Ao (long-dashed lines) all have shorter 
ages, as in the i — 1 case. The difference is that they are now even shorter, because the negative cosmological term 
is driven to high negative values more quickly in the past direction when £ = 2. Thus, the age of the Ao = — 1 model 
has dropped from 0.48 to 0.33 Hubble times, while that of the Ao = —3 model is now only ro — 0.30, down from 0.35 
in Fig. m. 

Fig. p[ indicates that the Ao < models tend toward eventual recoUapse, as they did in the £ — 1 case. However, 
this process now takes much longer. In other words, while the £ = 2 models are younger, their life expectancies are 
considerably greater. This can be understood by means of the same argument as before (§ 111 C| ). The larger value 
of £ means that the contribution of the cosmological term to the deceleration drops off more quickly in the future 
direction, thereby postponing recollapse for a longer period of time. Thus, models with Ao = —1 and Ao = —3 now 
survive for t S> 5 and t — 3.37 Hubble times respectively. 



E. The Case £ = 3 



For this case, eq. (3.7) reads: 



r —5- - ax = 0, 
dr-^ 



where a is given by eq. (3.10): 

a = (37/2)2Aor3. 

Following standard techniques [^9|, we switch independent variables from r to z = 1/r, whereupon: 



(3.41) 



(3.42) 
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d'^x ^dx 
dz^ dz 



ax = 0. 



(3.43) 



This is now in a similar form to eq. ( |3.11 ) in the I — 1 c ase, a nd it can be solved in the same manner Jzof . Changing 
independent variables again, from z io y = 2^— az, eq. (3.43) takes the form: 



d^x 
dy"^ 



dx 2 
dy 



0. 



(3.44) 



This is again transformable to Bessel's eq uati on, but with a different general solution x{y) — y ^[ciJi{y) + C2Yi(?/)], 
where ci,C2 are arbitrary constants. Eq. (|3.9D then gives a(r), as usual: 



a(r) 



- T-l/37 



[cMy) + C2Y,iy)f'\ 



(3.45) 



where y(r) = 2-\/— a/r and we have absorbed a factor of 2^— a into ci, C2. The Hubble parameter is found as usual 
by putting x(y) into eq. (3.6): 



li{T) - ifo 



c\Ji{y) + 02^2 (y) 



ci-^i(2/) +C2yi(y) 



(3.46) 



where 3i(y) and ^2(2/) are Bessel and Neumann functions of order two. As in the £ = 1 case, these solutions are 

real-valued (for t > 0) only if the cosmological t erm i s nega tive. 

In conjimction with the boundary conditions ( ^.16| ), eqs. (3.45) and (3.46) give for c\ and C2: 



ci 



C2 



V^i"2(yo)-ll(yo) 



V-Aqto [Ji (2/0)^^2(2/0) - -^2(^0)^1(2/0)] 
-V-Ao ^2(2/0) + >/i(2/o) 



V-AoTo [Ji (2/0)^1^(2/0) - -^2(2/0)5^1(2/0)]' 



(3.47) 
(3.48) 



where j/g = ^0 is given by eq. (3.19). We require one additional boundary conditio n to f ix tq. Unfortunately, as in the 
previous subcase, the procedure is complicated by the fact that both terms in eq. (3.45) diverge at t = 0, whereas we 
expect that the scale factor as a whole should behave smoothly there. 

We can make this more precise by employing the asymptotic expressions for Ji(2/) and Y\{y) at large y (ie., small 
t). We find (for t < 1): 



air) 



-I/27 



Ci sin 



^0 



C2 COS 



Wo_ 



^ 2/37 



(3.49) 



where C\ = Cq^^'^{ci 
for 1 

t ~2 case. Once tq is obtained, the values of ci and C2 follow from eqs. ( 3.47 ) and ( 3.48 ) respectively. 

The evolution of the scale factor for this case is illustrated in Fig. ^, which has exactly the same format as Fig. 
The Ao = model (solid line) has an age of 2/3 Hubble times, while those with negative values of Aq (long-dashed 
lines) all have shorter ages, as usual. The age of the Aq = —3 model, for example, has dropped from 0.30 to just 
0.27 Hubble times. It is u nlikely that any of these models could describe the real universe, given the observational 
constraints on Ho and To (j jlll C ). 

The main difference between Fig. |^ and its predecessors occurs at large times, where we observe that the curves 
all straighten out, and show no sign of leading to a recoUapse of the scale factor. The explanation for this is that the 
(negative) cosmological term is now decaying so quickly with time that, like ordinary matter, it is no longer sufficient 
to turn the expansion around. 



C2), C2 = —Cq^^'^{ci -f C2), Co = I/y^ttZJo and luq = 37x0%/— Aoto. This goes smoothly to zero 
0. In order to find the correct (ie., self-consistent) value of To, we use the sa me n umerical approach as in the 



F. The Case I 



For this case, eq. ( p. 71) reads: 



d^x 
dr^ 



^ a X ~ 0, 



(3.50) 
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with a given by eq. (3.10) as follows: 



(3.51) 



Switching independent variables via z = l/r as in the previous case, we find that eq. (3.50) takes the form: 

d^a; ^ dx ^ 
z—^ + 2—-azx^0. (3.52) 
az^ dz 

Employing standard methods ||79|], we switch dependent variables from z to y via dy/dz = zx{z). It may then be 
verified that y{z) satisfies the familiar equation: 



d^ 

dz^ 



ay. 



(3.53) 



The solu tions of eq. (3.53) are well-known; there are three subcases to consider, according as a [and hence Aq, 
eq. ( |3.5l| )] is positive, zero, or negative. 



1. The Subcase Aq > 



Since we exp ect o n observational grounds that Aq is probably positive, this is the most physical of the three. 
Solution of eq. ([3.53| ) for y{z) and hence x{z) is straightforward. The scale factor and Hubble parameter are given by 
eqs. (E|) and (pF~ 



air) 
H{t) 



r2/37 



2go 

37 



ci exp 1^ j + C2 exp 

ci(l - \/a/T) exp(\/a/r 



2/37 



C2(l + Va/r) exp(-\/a/r 



cir exp(-ya/T) + C2T exp(— ^/a/r) 

unless ci 



where ci,C2 a re arb itrar y con stants. It is clear that a(r) diverges at r 
assume. Eqs. and ( |3^ ) simplify: 



a{r) 
H{r) 



C2T exp 



2/37 



/a 



(3.54) 
(3.55) 

0, which we consequently 

(3.56) 
(3.57) 



37 V 

We then apply the boundary conditions ( |3.16 ) a,t t — tq, as usual. Eq. ( 3.57 ) with H{tq) — fixes the age at 

2 



To 



37(1 - VM) ■ 



(3.58) 



Substitu ting t his result into eq. ( 3.56 ) with a(To) = 1, we find that C2 — (I/tq) exp(-\/a/To), which can be put back 
into eq. ( 3.56 ) to yield this expression for a(r): 



a(T) 



— ) exp 

To 



a 
To 



-,2/37 



(3.59) 



We can draw a number of conclusions from eq. ( 3.58 ): firstly, that Ao < 1; and secondly, that the age of th e Uni verse 
To ^ 00 as Aq — > 1. This is reminiscent of the i = 2 case, and in fact eq. (3.58) is almost identical to eq. (3.27), the 
only difference being that Aq in the denominator has been replaced by \/Ao. Therefore, for the same value of Ao, the 
^ = 4 models are longer-lived by a factor of (1 — Ao)/(l — \/Ao)- This is due to the fact that, for positive A, the higher 
value of £ means that the cosmological term exerts a greater repulsive force in the past direction, and is consequently 
able to push the big bang back more effectively. As in the t — 2 case, we also find a lower limit on the age of these 
models. This is, from eq. (3.58|): 



To > 2/37. 



(3.60) 



This is twice as long as in the 1 = 2 case, but here the reason is simply that this lower limit corresponds to the case 
Ao = (not Ao = — 1 as before). With dust-like conditions (7 = 1) and the upper bound Hq < 83 km s^^ Mpc^^ [ pl| , 
we now find a minimum age of at least 7.9 billion years. 
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2. The Subcase Aq = 



For this subcase eq. ( 3.53| ) is trivial. Using eqs. (|3 

a(r) = 
H{r) = 



(c3^ 

2Ho 



I) and (|3 



we find immediately that: 



C4 



37 \C3 + C4T 



(3.61) 
(3.62) 



where 03,04 are arbitrary constants. For the first time we have a scale factor with the potential to go smoothly to 
some finite value other than zero at r = 0. Let us pursue this possibility and see if a nonsingular solution is possible. 
Instead of the boundary condition aijo) = 1, we impose a(0) = a*, where a* is the minimum value of the scale factor. 
In eq. (3.61) this implies that C3 = a^''^'^. Inserting this result back into eq. ( [3.61 ) and applying the usual condition 
a(To) — 1, we obtain C4 = (1 — C3)/to. Substituting this into eq. ( 3.62 ) and applying the third boundary condition 
H{to) = Hq, we obtain for the age of the Universe: 



To - 2(1 - C3)/37. 



(3.63) 



This result matches onto that of the previous subcase, eq. (3.58), only if 03 = 0, which also implies that a* = 0. The 
present case is therefore singular at r = 0, like all the others studied in this section. The age of the Universe is given 
by eq. ( 3.63| ) as tq = 2/37. Putting these results back into eq. (B.61), we find: 

2/37 



air) 



To 



(3.64) 



This is just the standard fc = solution with no cosmological term, as might have been expected. 



3. The Subcase Aq < 



Solution of eq. ( |3.53| ) is again straightforward and leads via eqs. ( |3.9D and ( |3.6| ) to: 

2/37 

C5 sm I I + cg cos f 



a(r) = r2/3^ 
iJ(r) 



T / V T 

2ffo f [c5 + (\/^/r)c6] sin(V^/T) + [ce - {^f^/T)o5] cos(v^/r) 



(3.65) 



37 



05T sm{^/— a/ t) + cqt cos(-\/— a/r) 



where 0^,0^ are arbitrary constants. Application of the boundary conditions (3.16) fixes these constants in terms of 
To- 



1 



C5 



C6 



37 

Po sin ^0 + cos ^0 ( 1 - Y 
37 

Po cos (3o - sin ^0 ( 1 - 



(3.66) 
(3.67) 



where /3o = (^) toV~^()- As usual, we require one additional boundary condition to fix the value of tq. The situation 
is again complicated by the fa ct that both terms in eq. (3.65) diverge at r = 0, whereas ajr) itself goes smoothly to 
[In fact, eq. (3.65) has exactly the same form as the asymptotic expression ( 3.49| ) in the £ = 3 case.] We 

Once Tq is obtained in this way, the 



[II D 3 



zero there. [In fact, eq. (3.65) has exactly the same form as the asymptotic expression (3. 
therefore have recourse once again to the n umeri cal m ethod described in 
values of 05 and cg are fixed by eqs. ( |3.66 ) and ( 3.67 ). 

The evolution of the scale factor for this case is illustrated in Fig. which has exactly the same format as Fig. |. 
Several features may be noted. To begin with, we see that models with A cx are qualitatively the same as those 
with A oc for positive Aq, and qualitatively similar to those with A oc t^^ for negative Ao. 

There are important quantitative differences, however. Models with positive Aq are significantly older. With 
Ao = 0.5, for example, tq is now 2.28 Hubble times — older by a factor of 1.71 times than the equivalent £ = 2 model, 
exactly as predicted in the discussion following eq. (3.59). This is once again due to the fact that, with ^ = 4, a 
positive A-term increases in size very rapidly in the past direction. Negative- Aq models, on the other hand, are once 
again younger. The age of the Ao = —3 model, for instance, has dropped from 0.27 to just 0.25 Hubble times. And 
finally, in the future direction, we see that there is no longer very much distinction between the Aq < and Aq = 
models, compared to Fig. ||. The contribution of the cosmological term to the deceleration of the scale factor now 
drops off so quickly that it rapidly becomes irrelevant. 
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IV. A AS A FUNCTION OF THE SCALE FACTOR 



A. Previous Work 



We now move on to consider decay laws of the form set out in eq. (2.11): 

A = 6a-'". (4.1) 

The scale factor may be more natural than time as an independent variable, to the extent that many physical quantities 
(such as temperature) depend more simply on a than t. Nearly half of the decay laws listed in Table | contain terms 
of the form ( |4.1| ) . The case m = 2 has been singled out for the most attention | |4l| , |4^ -^ , and is motivated by some 
dimensional arguments |^,^. A second group of authors has concentrated on values of to « 4 [^-^, for which 
the A-term behaves like ordinary radiation. It has been shown that, for certain kinds of A-decay, the lower of these 
two values of m is thermodynamically more stable |^^. A third idea, that the cosmological term scales with a like 
ordinary matter (to = 3), follows from one interpretation of an intriguing new scale- invariant extension of general 
relativity [Q. 

There are also some studies in which the value of m is not fixed a priori. Ages of general-m models have been 
calculated and agree with observation if to < 3 . The power spectrum of matter density perturbations does not 
appear to be greatly modified by a decaying A-term, at least for < to < 2 . Lensing statistics favor models with 
TO > 1 1^^, or TO > 1.6 when combined with other tests involving CMB anisotropics and the classical magnitude- 
redshift relation for high-rcdshift supernovae Other aspects of models in which A decays as a^™ have been 

discussed by several authors p3| , ^j59| , although no specific numerical bounds are set on m. 

The question of the initial singularity has so far received little attention in theories of this kind . Some explicitly 
nonsingular solutions have been constructed, all with m = 2 |Q. In one other case it is noted in passing that the 
existence of an initial singularity would require < m < 4 under some circumstances pof . The remaining authors 
either do not mention the issue, or (as in one case ||5l|| ) rule out a priori the possibility of nonsingular solutions. In 
this paper, we take a broader view and examine all possible solutions for the scale factor, including those in which it 
takes a nonzero minimum value. Moreover we will extend the discussion, not only to general m, but to general 7 as 
well [where 7 characterizes the equation of state of ordinary matter, eq. ( ^.6D ]. If to and 7 are thought of as defining a 
parameter space, then we determine, firstly, the extent to which the space is singularity- free; and secondly, the extent 
to which it is observationally viable. 



B. Evolution of the Scale Factor 



We begin with the dynamical equations (2.4) - (2.9). In particular, we consider the energy conservation law (p. 7]), 
which, with the decay law (11), becomes: 



37^ _ 



8^ 



,37-(™-|-l) 



Integrating, we find for the matter energy density: 

p{a) = poa^^'^fia), 
where we have set a{to) — ao — 1, p{ao) — po, and defined: 



/(a) = 1 + Ko X j 
kq = B/SnGpo- 



37 — m 
37 ln(a) 



if TO 7^ 37 
if TO = 37 



(4.2) 
(4.3) 

(4.4) 
(4.5) 



When TO = 0, then /(a) = 1 and eq. (12) yields the usual result that p scales as in a pressure-free universe (7 = 1) 
and a^* in a radiation-dominated uni vers e (7 = 4/3). The new parameter kq can be fixed in terms of observable 
qua ntitie s by means of the decay law (11), which gives S = Aq = '3H§Xo [with uq — 1 and Aq defined as usual by 
eq. ( ^.15| )]. Substituting this result into eq. ( f4.5| ), we find: 



Ko — Ao/f^Oj 



(4.6) 
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where is defined by eq. ( ^.15| ) . The parameter kq is simply the ratio of energy density in the cosmological term to 
that in ordinary matter at the present epoch. 



Substitution of eqs. (1.1) and (4.3) into the Lemaitre equation (2.4) yields 



da 
d^ 



a [Qoa^^^'fia) - i^o + Aq - l)a~^ + Xoa~ 



I 1/2 



(4.7) 



where we have made use of the definitions ( ^.15| ), recaUed that a/a = {Ho/a)da/dT, and selected the positive root 
since redshifts rather than blueshifts are observed. 



At this point we could choose integer values of m and attempt to solve analytically for a(T), as in §111. Detailed 
analyses have been carried out along these lines for the case m — 0] ie., for a constant cosmological term |8^. It is 
doubtful that they can be usefully extended to the general situation in which m 7^ 0. We opt instead to solve the 
problem numerically, following the lead of Felt en and Isaacman ||^. The time-derivative of eq. ( [1.7| ) is: 



d^a 



noa'-'yia) + \oa'- 



(4.8) 



[This could equally well have been obtained from eq. (| 
for the scale factor: 



8h.] We substitute eqs. (4.7) and (4.8) into a Taylor expansion 



da\ ^ 



(4.9) 



This can be integrated numerically backwards in time to determine whether or not a model with given values of 
{to, 7, fig, Ao} eventually reaches a — 0. We have tested the procedure for the case of a constant cosmological term 
(to — 0) and dust-like equation of state (7 = 1), and our results in this case confirm those of Felten and Isaacman [|]. 
Fig. H depicts a group of examples with ilp — 0.34 (typical of large-scale observations |^^). and various values of Aq 
(labelled beside each curve). Note that the difference between this figure and the ones in the previous section is that 
we now include models of all three kinds: closed (long-dashed lines), flat (dash-dotted lines) and open (short dashes). 
(To keep the diagram from being too crowded, we show only models with the same value fig of the matter density.) 
Fig. ^ indicates that negative values of Aq can lead to recollapse in open, as well as flat universes (c/. §|l|). Of greater 
interest, however, is the fact that models with Aq above a critical value A* (=1.774 605 in this case) avoid the big 
bang, undergoing a finite "big bounce" instead. Models with slightly less than this critical value (eg., Aq = 1.76 in 
this case) are of the "coasting Lemaitre" kind: they begin in a singular state but go through an extended phase in 
which the scale factor is nearly constant. Models with exactly Aq = A* (shown in Fig. |^ with a solid line) are perhaps 
the most interesting of all. As time r — s- —00, they neither plunge to zero size nor bounce back up to infinite size, 
but level off indefinitely at a constant value a = a* (=0.46 in this case). These are nonsingular Eddington-Lemaitre 
models, asymptotic to the static Einstein universe in the infinite past. 

All these features of the m = 0, 7 = 1 models have been discussed at greater length elsewhere [||. Our purpose here 
is to generalize the discussion to arbitrary values of m and 7. 



C. Critical Values of Ao 



In particular, we wish to obtain general expressions for the critical value A* of the lambda parameter Aq and 
minimum size a* of the scale factor, given any class of models {m,7, fio}- As discussed in §[IIA, we are interested 
in models for which da/ dr ^ at some time in the past. This occurs, for example, at the moment of the "bounce" 
in all the oscillating models shown in Fig. ^. The critical case is distinguished by that fa ct th at no t on ly da/dr, but 
also d'^a/dr^ vanishes at this point We therefore set da/dr = d'^a/dr'^ = in eqs. (^) and ( |4.8| ). This yields: 



A* 



(37-2)(37-TO)f}o 



37(2 - TO)af + (37 - 2) 



TO 



(4.10) 



where a,, the minimum value of the scale factor, is found by solving: 

37(37 - TO)r!oa^ -I- 37(2 - to)(1 - no)al-< - 

-(37 -2) (37^*0 - m)ar = 0. (4.11) 

In general eq. (4.11) has to be solved numerically, but in the case to = 0,7 = 1 it reduces to a cubic equation 
Fig. g is a phase space portrait of this case, with each point on the diagram corresponding to a choice of VIq and 
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Ao (after Lahav et al [|T^). The critical values A* are represented in this figure by a heavy solid line. The region 
to the right of this curve corresponds to universes with Aq > A*; that is, with no big bang. Also shown in Fig. g is 
a straight dash-dotted line representing the boundary between open and closed universes ; mo dels on this line have 
k = (r^o = 1 — Ao) while those on the left and right have fc < and fc > respectively (§ [ID ). 

We now have the tools we need to investigate models with arbitrary values of m and 7. The idea will be to use 
phase space diagrams like Fig. ^ to determine how much of the parameter space (1) corresponds to models with a 
nonzero minimum scale factor; and (2) agrees with observational constraints. We can then confirm whether a given 



model with {m,7, Qq} does in fact avoid the big bang by carrying out the numerical integration described in §[VB, 
and plotting the results on evolution diagrams hke Fig. |. 



D. Observational Constraints 



1. Upper Bounds on Ao 



We pause first to take stock of some of the experimental constraints that have been placed on models with nonzero 
cosmological terms. Most immediate are direct upper bounds on Aq from a variety of methods, most of them assuming 
that JIq + Ao = 1. Until recently, these have typically been of order ~ 1 |]l5|. Additional methods, however, have 
become available in the past few years. CMB fluctuations, for instance, have produced an upper limit of Ao < 0.86 [Q. 
Gravitational lens statistics give a tighter bound of Ao < 0.66 and observations of Type la supernovae appear 
to reduce this still further, to Ao < 0.51 [Q. All of the above are described as 95% confidence level measurements. 
On the other hand, a lower limit of Ao > 0.53 has been obtained from the galactic luminosity density — also at 95% 
confidence po|| . 

Complicating the picture somewhat is the fact that several other observational data are well-explained by substantial 
values of Ao, including the lack of observed small-scale dark matter, the expectation that infiation should lead to 
near-fiatness, and especially the high age of the Universe inferred from models of stellar evolution. This "cosmic 
concordance" at one time led to calls for Ao '■^ 0.8 

Thus, the observational situation is not yet settled. It is clear, however, that the nonsingular models in Fig. |^, 
which require Ao > 1.77, are almost certainly unphysical. It remains to be seen if the same conclusion applies when 
TO ^ and/or 7 7^ 1. 



91|, although this has since dropped to 0.5 — 0.7 195 



2. Age of the Universe 



A lower limit on Aq derives from the age of the Universe, — Jq° da/a. In our units of Hubble times: 

da 



To = Hoto 



da/dr 



(4.12) 



If we use a recent value of 73 ± 10 km s~^ Mpc^^ for Hq in conjunction with the oldest globular cluster 
age of to > 9.6 Gyr |^^, then tq > 0.62; that is, the Universe is at least 0.62 Hubble times old. Of course, as 
noted in § 111 C| , the true status of these two parameters is still a subject of some controversy. A lower value of 
Hq 55 ± 10 km Mpc~^ [|3| leads to tq > 0.44. We choose an the intermediate value of 0.5. In conjunction with 
our expression (4.7) for da/dr, this implies: 



To 



dv 



> 0.5, 



'1 [nov^+^ygiv) - {rto + Ao - 1)^4 + Aoa2+™]^/2 
where we have changed variables of integration from a to v = a^^ = 1 + z for computational reasons, and: 



(4.13) 



g{v) = 1 + Ko X \ 37 — TO 
{ — 37ln(?;) 



if TO 7^ 37 
if TO = 37. 



(4.14) 



Eq. (4.13) reduces to the sta ndard results in the case to = 0,7 = 1. 

Numerical solution of eq. ( 4.13 ) produces a lower limit on Ao as a function of {to, 7,^0}- This age constraint is 
shown in Fi g, (for the case to = 0,7 = 1) as a long-dashed line. Its position matches that in a similar plot by 
Lahav et al [psf. The region to the left of this curve corresponds to universes younger than half a Hubble time. Big 
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bangless models, of course, are not constrained by this; they are infinitely old (by definition), the heavy solid line 
being precisely the boundary where tq — > oo. The main impact of the age constraint is to rule out models with a 
negative cosmological term. 

In Fig. 0, we show the effects of varying the parameters m and 7 respectively on this age constraint. It may be seen 
that altering the value of m changes the slope of the curve, but does not otherwise greatly affect the age, even over 
the range — 1 < to < 3. Altering the value of 7, on the other hand, has a larger effect. In particular, the "harder" 
the equation of state (ie., the larger the value of 7), the further this constraint encroaches on the available parameter 
space. This is in accord with the well-known fact that a radiation-dominated universe (7 = 4/3), for example, is a 
short-lived one. 



3. Gravitational Lensing and the Antipode 



For closed models, the most stringent constraint on Aq comes from gravitational lensing, which requires that the 
"antipode" be further away than the most distant normally lensed object The antipode is the point where the 
radial co ordin ate u) = dt/a{t) = da/aa takes the value tt Using our expression (4.7) for da/dr, together 
with eq. (2.16), we obtain (taking ao = 1 as usual): 



v/f^o + Ao - 1 



dv 



[Qov^^giv) - {Qo + Ao - l)^^ + Aoa" 



11/2' 



(4.15) 



where we have again changed integration variables from a to v. The standard formula g] is recovered when m = 
and 7 = 1. 

At present, the furthest known normally lensed object is a pair of lensed galaxies at Zi — 4.92 ||9^. We therefore 
require: 



ZA(f^o,Ao,TO,7) > 4.92, 



(4.16) 



where za{^o, Xo,m,'j) is defined by eq. ( 4.15| ) with w = tt. Numerical solution of this equation yields an upper 
limit on Aq as a function of {m,j,QQ}. This lensing constraint is shown in Fig. |^ (for the case to = 0,7 = 1) as 
a short-dashed line. Its position is close to that in previous plots [^,^ employing a smaller value ze = 3.27 (our 
constraints are slightly stronger). The region to the right of this curve corresponds to universes incompatible with 
the lensing observations. Since this includes the entirety of nonsingular parameter space, we can see that oscillating 
models with constant A and zero pressure are ruled out, as noted previously by these authors. 

In Fig. 1^, we show the effects of varying the parameters m and 7 respectively on this lensing constraint. As before, 
it is seen that harder values of 7 lead to tighter constraints on the available parameter space. However, the situation 
with regard to to is altered quite dramatically. In particular, the higher the value of to, the weaker the lensing 
constraint becomes. As we will see, this significantly improves the prospects for viable big bangless models. 



4- The Maximum Redshift Constraint 

A fourth observational constraint, which must be satisfied only by nonsingular models, concerns the maximum 
observable redshift = — 1 in a universe with a minimum scale factor a,. This must obviously be larger than 
the greatest redshift Zobs actually observed. Thus the nonsingular models in Fig. ^, which never get smaller than 
a^, = 0.46, cannot accommodate redshifts greater than z» — 1.2. This disagrees with observations such as those of the 
lensed galaxies mentioned above. 

The condition > Zots can be reformulated as an upper limit on the m atter d ensi ty of the Universe. Assuming 



that da/dr = and d^a/dr^ > at z = z, for nonsingular models, eqs. (4.7) and (4.8) can be combined to read (for 
TO 7^ 37): 

^ 2-TO+(TO/37)(37-2)(l + z,)37-r" 

° - 2 - TO - (37 - to)(1 + z,)3-/-2 + (37 2)(1 + z,)3t-™ ■ ^ ^ ^ 



Eq. (4.17) reduces to earlier expressions of Borner and Ehlers ]l4| when to = and 7 = 1 or 4/3. These authors have 



then argued as follows: given that quasar redshifts have been observed out to Zohs > 4, we know that the Universe 



has z* > 4. This constraint with eq. (4.17) implies (assuming m = and 7 = 1 or 4/3) that flo < 0.018, which is 



contrary to observation. This indicates that our Universe could not have been nonsingular. 
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Let us see how the above conchision changes when we generahze the situation to values of m ^ 0. Moreover, 
we will strengthen the argument by noting that some distant galaxies have now been assigned photometric redshifts 
as high as Zobs > 6 implying that > 6. The resulting upper limits on density are listed in Table || for 
various values of m and 7. From this table we see that the new photometric redshifts tighten the Borner-Ehlcrs 
constraint noticeably: as long as m = 0, a nonsingular Universe requires flo < 0.006 (if 7 = 1) or < 0.001 (if 
7 = 4/3). These numbers, of course, are too low to describe the real universe. However, Table || demonstrates that 
much higher densities are possible in singularity-free models with variable A. For example, retaining pressure-free 
conditions (7 — 1), we see that if to = 1, then the matter density in a nonsingular universe must satisfy f2o < 0.48. 
This value is not unreasonable at all; in fact it is well above most dynamical measurements, which suggest J7o ~ 0.3 
p7| . The constraint is similarly loosened if we move toward softer equations of state with 7 < 1, such as those that 
have been proposed in |3^ . 



E. Viable Oscillating Models 

We now demonstrate that models with to ^ are capable of satisfying all the constraints discussed above. Fig. ^ 
is a phase space diagram like Fig. |[ but plotted for a range of nonzero values of m. We have assumed 7 = 1 as 
before, but there is now a different line of critical values A*(m,7, fio) for each value of to. Like Fig. |^, Fig. || shows 
that much of the nonsingular parameter space (below and to the right of the heavy solid lines) is eliminated because 
it does not overlap with the regions allowed by the lensing constraint (above and to the left of the lighter dashed 
lines). This is especially true for small values of to. With increasing to, however, significant triangle-shaped regions 
of parameter space appear near fc = 0. Thus for to = 1, there are allowed models with Qq ^ 0.3 and Aq between 
about 0.7 and 0.9. These are very close to the values favored by observation so we focus on this case. The 

lensing constraint imposes the lower bound fio > 0.31. If we take flo = 0.34 as a specific example (as before), then 
by tracing horiz ontall y across Fig. ^ we find that lensing also places an upper bound Aq < 0.72 on the cosmological 
term. From eq. ( 4.10|) the critical value for this case turns out to be A* = 0.68, which is marginally consistent with 



the upper limits on Aq mentioned earlier. Therefore to = 1 models with Oq = 0.34 and 0.68 < Aq < 0.72 are both 
realistic and singularity- free. 

For larger values of to, the range of acceptable Ao-valucs is broader, but one is also driven to higher values of fig. 
With m = 3/2, for example, we find that viable nonsingular models occur only for fig > 0.46, and that at JIq = 0.51 
they lie in the range 0.51 < Ao < 0.57. At m = 7/4 and flo = 0.59 this broadens to 0.42 < Ao < 0.50 (low enough to 
satisfy even the supernova constraint |^^ ). As to — > 2, we find that a* — > 0, so that nonsingularity is lost. This also 
happens in the limit 7 — > 2/3 and, interestingly, as A* ^ 1 — JIq. (It is for this latter reason that all the nonsingular 
solutions in Fig. ^ have k > 0.) 

Unusual values of 7, although interesting and possibly relevant at early times, do not alter these results in any 
fundamental way. In practice we find that the parameter space allowed by observation shrinks slightly for harder 
values of 7, like 4/3, but grows significantly for softer values such as those in the range 0.4 < 7 < 1 |p7|. 

To verify that models with the properties described above can in fact avoid the big bang. Fig. IG shows the 
evolution of the scale factor when to = 1,7 = 1 and fig = 0.34, obtained by numerical integration using the Taylor 



expansion (4.9), as before. Various values of Aq are labelled beside the appropriate curves. This plot has exactly 
the same format as Fig. |^, with long dashes corresponding to fc > 0, short ones to k < 0, and the dash-dotted line 
corresponding to fc = 0. 

Fig. ^ confirms that when Aq takes on the critical value A* (=0.679 807 621 in this case), the scale factor a 
evolves back to a constant value a*, as expected (solid line). More importantly, the small size of this minimum value, 
a* = 0.0097, means that we can now accommodate observed redshifts up to = 102, well beyond the furthest objects 
yet seen. The model is thus compatible with all observational data. It cannot, however, accommodate larger redshifts 
like that attributed to the last scattering surface {ziss ^ 1100), let alone the era of nucleosynthesis, which occurred 
at temperatures Tnuc ^ 10^° K, or — since T oc (1 -I- z) — redshifts Znuc ^ 10^°. In fact, with z^, = 102 the m — 1 
model heats up to no more than about = 103ro = 281 K, which is just under room temperature. That is to say, 
if the Universe is described by a nonsingular variable- A theory with to = 1 , then the CMB and the abundance of the 
light elements cannot be accounted for in the conventional way. 

Unconventional explanations have been proposed for both these phenomena, but are not widely accepted [p6| . A 
more conservative approach might be to retain the traditional hot, dense early phase in the context of variable- A 
theory by moving to larger values of to. As discussed above, this leads to smaller values of a* and hence larger 
values of z* and T*. In this picture the high temperature of the early universe is a result, not of an initial big bang 
singularity, but of a very deep "big bounce." With m = 1.5 and JIq = 0.51, for instance, we find that a* = 1.7 x 10~*, 
corresponding to a maximum redshift of 5900 and a "bounce temperature" of 16,000 K — more than sufhcient for 
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recombination. At to = 1.8, this latter number climbs to 5.7 x 10^ K — hot enough for nucleosynthesis. (An analysis 
like that below Fig. ^ shows that models with to = 1.8 and fig = 0.61, for example, are observationally viable if 
0.41 < Ao < 0.49.) With to — 1.9, one can entertain bounce temperatures as high as 2.9 x 10^^ K (2.5 GeV), 
approaching the realm of elcctroweak unification. (Models with to = 1.9 and Qq — 0.67 agree with observation if 
0.39 < Ao < 0.48.) As TO 2, in fact, one finds that a* — *■ and oo. As remarked above, however, these higher 

values of to come at the modest observational price of higher matter densities. The lensing constraint sets a limit of 
fio > 0.54 for TO = 1.8 models, for example, and ilo > 0.57 for to = 1.9 models. 



V. A AS A FUNCTION OF THE HUBBLE PARAMETER 



A. Previous Work 



appear to offer the possibility of avoiding the big 
The parameter space occupied by viable models, 
one reason for this is that the oscillating solutions are all closed, 
whereas observational evidence tends to favor an open universe |p6| . In this section we shift our attention to decay 
laws of the form set out in eq. (2.12): 



We have seen in §IV that variable-A theories in which A cx a^" 
bang singularity without violating any obs ervat ional constraints, 
however, remains small. As discussed in SIVE, 



A = CH' 



(5.1) 



Among other things we will find that a cosmological term of this kind allows for a nonzero minimum scale factor in 
open models. 



Terms of the form (5.1) are nearly as common as those studied in ^V, appearing in almost a third of the decay laws 
listed in Table |. The value n = 2, favored on dimensional |Q and other grounds |^^, is adopted almost universally 
Other powers of H are considered in three cases, but only in combinations such as exp(— J aH"' dt) 
[70[, and aH{dH/da) The question of the initial singularity has received little attention in these 

papers, apart from one explicitly nonsingular and one asymptotically de Sitter-like solution ]70|]. 




B. Riccati's Equation 



Substituting the decay law (5.1) into the differential equation (|2.9|), we obtain: 



dH 
da 



37-2 
2a 



(5.2) 



where we have us ed the fact that d/a — aH{dH/da) + H^. This is nonlinear, but bears some resemblance to the 
Riccati equation ( |3.5| ). We switch dependent variabl es fr om H to x = , where r is an arbitrary constant whose 
value will be chosen in a moment. In terms of x, eq. ( |5.2| ) takes the form: 



dx / r^C\ ( / 3r7 



da V 2a 



V 2a 



r(37- 2)fc' 



2a3 



where £ = {n + r — 2)/r and to = (r — 2) /r. 

When n = 4, we see that the choice r = 2 puts i — 2 and to = 0, whereupon eq. (^) takes the form: 

— = V{a) x^ + Q(a) x + 7^(a), 
da 



with: 



V{a) = — , Q(a) = , 7^(a) = ^ . 



(5.3) 



(5.4) 



(5.5) 



This is Riccati's equation (3^), as desired. 

When n = 2, on the other hand, the choice r = 2 leads io I = \ and m = 0, whereupon eq. ( ^.3D reduces to: 

^ + 5(a) X = T{a), 



(5.6) 
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with: 



_ 7(3 -C) 



r(a) 



_ (2 - 37)fc 



(5.7) 



This case is hnear, and may be solved easily. 

For other values of n, a check reveals that no apparent choice of r will put eq. (x3) into a form which can be solved 
analytically for x (ie., for H). We therefore concentrate on the cases n = 2 and n = 4 for the time being, leaving the 
others to future numerical analysis. In both cases the Hubble parameter is given hy H — y/x (since r = 2). 



Multiplying through by a factor of exp [/ S{a)da\ = a 
directly for x{a) and hence H{a): 



C. The Case n = 2 

7(3-c) p^^g gq^ (5.6) into exact form, which may be integrated 



H{a) 



(2-37)fc 



1/2 



(5.8) 



_7(3-C)-2j 

where Cq is a constant of integration. Using the boundary condition H^uq) = H[l) = Hq to eliminate Cq, we obtain: 

1/2 



where: 



H{a) = Ho 



,7(C-3) 



Pa' 



2-37 



7(3 -C)- 2 



Hi 



(5.9) 



(5.10) 



The parameters a, (3 and C can all be fixed in terms of observable quantities at the present time t = to, as follows. 
With n = 2, the decay law (|5.lD gives C = Aq/Hq = 3Ao, where we have used the definition (2.15). Substituting this 
result into eqs. ( ^.10| ) and (^.9|), we find: 



a 



(37/2 - i)no - Ao 



(37/2)(l-Ao)-l 
da 



dT 



a a 



2-37(l-Ao) 



, f3=l-a, 
1/2 



(5.11) 
(5.12) 



where we have used eq. ( 2.16 ) and set oq = 1 as usual. 



D. Evolution of the Scale Factor 



For special values of 7 and Aq, it may be possible to solve eq. ( ^.12| ) in terms of elliptic (or simpler) integrals. For 
general purposes, however, we take the same approach as in §[V and evolve the scale factor numerically in terms of 
its first and second derivatives. The latter of these is given by: 



fa 



1-|(1-A„) 



l-37(l-Ao) 



(5.13) 



Eqs. ( p^ and ( jsl^ ) can be substituted into the Taylor expansion ( |4.9| ) and integrated backward numerically to 
determine the behavior of the scale factor. 

We also wish to determine the conditions under which the Universe evolves backward to a nonzero minimum scale 



factor, a — a*. Setting da/dr = (at a — a*), we find from eq. (5.12) that: 



(5.14) 



It is conv enient to distinguish two cases, according to whether the exponent on a* is positive or negative. Using 
eqs. (5.11), we obtain: 
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2Ao - (37 - 2)»o 



(37- 


2)(1- 


^0 — Ao) 


(37- 


2)(1- 


^0 — Ao) 


2Ao 


-(37 


-2)r!o 



l/[37(l-Ao)-2] 



l/[2-37(l-Ao)] 



if Ao < Ac 



if Ao > Ac 



(5.15) 



where we have defined Ac = 1 — 2/37 (leaving the case Ao — Ac aside for the time being). 

From eq. ( ^.15 ) we draw a number of important conclusions: (1) Spatially flat solutions {Qq + Aq = 1) have either 
a* = 00 (if Ao < Ac) or a, = (if Ao > Ac). The former case is not interesting. The latter case is de Sitter-li ke, wit h 
the initial singularity pushed back into the infinite past. We have encountered this kind of solution before (§ [IID3| ). 

(2) With the modest assumption that 7 > 2/3 (ie., normal, non- inflationary matter), then we notice that closed 
solutions (ie., Ao > 1 — fio) must satisfy Ao < (37/2 — l)ilo, while open ones (Ag < 1 — fio) obey Ao > (37/2 — l)rio. 
These conclusions follow from requiring that the terms in square brackets be positive; ie., from the requirement that 
a* be a real number. (The exception in which the exponent is an even integer p occurs only for special values of the 
lambda parameter, Aq = 1 — 2/37 — l/6p7, and will not be considered further here.) 



(3) Requiring that < a, < 1, we learn that both the numerators and denominators in eq. ( 5.15 ) must be positive. 
Comparing their relative absolute magnitudes, we distinguish two possibilities: (a) if Aq < Ac, then both the numerator 
and denominator must be positive, since otherwise (37 — 2)(1 — fio — Aq) < 2Ao — (37 — 2)J7o, which reduces to Ao > Ac, 
contrary to the hypothesis. On the other hand, (b) if Ao > Ac, then we also find that both numerator and denominator 
must be positive, since otherwise (37 — 2)(1 — Qq — Ao) > 2Ao — (37 — 2)rio, which reduces to Ao < Ac, again contrary 
to the hypothesis. Therefore both the numerators and denominators are positive in all cases. 

(4) It follows from the conclusion (3), in conjunction with the assumption 7 > 2/3, that: 



An < 1 - On 



(5.16) 



In other words, to realistically describe the present Universe, models with a nonzero minimum scale factor must, in the 
present theory, be open. While some nonsingular open solutions have been found in theories employing scalar fields 
and higher-order curvature terms JOTj , we are not aware of precedents for this in theories based on the cosmological 
term. 



(5) By combining eq. (5.16) with the conclusion (2), we infer that: 

Ao > (37/2 - 1)^0, 

which sets a low er limi t on the size of the cosmological term. 

(6) Eqs. ( 5.16 ) and ( 5.17 ) together impose an upper limit: 



Oo < 2/37, 



(5.17) 



(5.18) 



on the matter density of the Universe. 



E. Minimum Values of the Scale Factor 



The information contained in eqs. (^.15 ) - ( 5.18| ) is summarized in Fig. (ITTl), which is a phase space plot like 



Fig. ||, but enlarged to show only the range of interest, < Oq < 1 and < Ao < 1. Models with (Ao,r2o) are 
represented by points on this diagram, as usual. The critical values of Aq in this theory define the upper edges of the 
triangular region at the base of the diagram; ie., the region bounded by the curves fio + Ao < 1 (dash-dotted Hne) 
and Oq + Ao > (37/2)r2o (dashed line). All models between these curves are nonsingular, with real values of a* in the 
range < a, <J_, as we have stipulated. Fig. |l^ is plotted for 7 = 1. 

Using eq. ( |5.15| ), we have plotted contours of equal minimum scale factor a* in this region (heavy solid lines). Any 
point along o ne of the contours corresponds to an oscillating model with the labelled value of a, . Following the 
discussion in § [V E| , we would like to find models in which a* is as small as possible, in order to obtain the largest 



possible maximum observable redshift Zobs < o,~^ ~ 1. For instance, to be compatible with observations of quasars 
{zobs ~ 10), a, must be less than about 0.1. If we wish to explain the CMB as relic radiation from the last scattering 
surface at zi^g « 1100, then we require a smaller minimum scale factor, a* < 0.001. And to meet the demand that 
the early Universe heat up to nucleosynthesis temperatures [T„„c ~ Tq (ao/a,) w 10^" K], our model must satisfy: 

a, < 10^^ (nucleosynthesis) . (5.19) 

Fig. |ri| demonstrates that the present theory can readily satisfy this constraint. Any model lying along the curve 
labelled 10~^ will be capable, in principle, of reaching these temperatures near the "big bounce." This confirms 



20 



comments made by several authors that there is no reason in principle why the oscillations in a nonsingular 
model cannot be deep enough to account for all the evidence which is usually taken as proof that the Universe began 
in a singularity. Moreover, the only observational constraint which seriously limited the variable-A models in the last 
section — the lensing cons train t — does not apply in this section because there is no antipode in an open universe. 
To meet the condition ( ^.19 ), models in this theory must lie close to the upper edge of the triangular region in 



Fig. 0. To an extent this is "fine-tuning." However, it also allows us to make very definite predictions (as in §rV) 
about the values of Aq that would be required in a realistic oscillating model. As an example, let us consider the 
observationally favored value of Hq — 0.3 and let us assume 7 = 1 as usu al. T racing horizontally across the line 
defined by flo = 0.3 in Fig. |ll], we can see that the nucleosynthesis condition ( |5.19| ) is met by only two values of the 



cosmological term: Aq « 0.15 and Aq « 0.7. Of these, the larger value is only marginally viable, being very close 



to the observational upper bounds described in §IVD1. The smaller value, however, is perfectly acceptable from an 
observational standpoint. 

In general, the theory predicts that the most likely value of Aq is either just below. 

1-no, (5.20) 

or else just above: 

(37/2-l)f^o; (5.21) 

which is to say, just above r2o/2 in a dust-like universe (7 = 1). The former situation might be preferable to some on 
theoretical gr ounds l|92| , while the latter is in better agreement with the increasingly stringent observational upper 
limits on Ao ( gIVDl| ). 

To confirm that models with these features really do avoid the big bang, th e so lutions can be evolved backward in 
time as before, using eqs. ( |5.12D and (|t|) with the Taylor series expansion (|4.9|). The results of this procedure are 
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shown in Fi g. ^2] for the case VLq — 0.3 and 7 = 1. This diagram is an enlarged version of the evolution plots in §IV 
(Figs. 1^ and |10|), showing only the past three Hubble times. Values of Aq are marked beside the appropriate curves. 
Several features can be noted. 

Firstly, the initial singularity is avoided for any value of Aq between 0.15 and 0.7, as expected on the basis of the 
phase space diagram. Fig. In the limiting case where Ao = 0.7 exactly, which is spatially flat, we see that a* = 



(de Sitter-like behavior), as expected based on the discussion following eq. ( 5.15 ). 

Secondly, Fig. |l2| confirms that the value of a* is smallest near the critical values of Aq: 0.064 for the Aq = 0.6 case 
(just below 0.7), and 0.018 for the Aq = 0.2 case (just above 0.15). These numbers have been chosen for illustrative 
purposes; smaller values of a* (with consequently larger bounce temperatures) are obtained by letting Aq approach 
the critical values more closely. 

Thirdly, this evolution plot gives us some information about the ages of the models; that is, the elapsed time since 
the big bang (or the big bounce, as appropriate). It may be seen that, within the range 0.15 < Ao < 0.7, larger values 
of Ao correspond to older universes, as usual: 1.10 Hubble times for Aq — 0.2, and 1.70 Hubble times for Ao = 0.6. 
Even if Hq takes on its largest currently acceptable value of 83 km s~^ Mpc^^ |^l|, the ages of the two models are 
13.0 and 20.1 billion years respectively — well above the globular cluster limit of 9.6 Gyr p2[. 



F. The Case n 



We proceed to the other case of interest, n — A, which consists of the Riccati equation (5.4) for x{a). This can be 
solved using standard techniques [[79[. Switching dependent variables from x to y via x — {—1/Vy) dy/da, we obtain: 



where: 



d^y . _idy _4 
-T^+C,a -r+m y = 0, 
aa^ da 



C=l + 37 , r7 = 7C(2-37)fc. 



(5.22) 



(5.23) 



The parameter ry can be connected to observation as fol lows. With n — 4, the decay law (5.1) gives C — Aq/Hq 
SXq/Hq. Substituting this result into the second of eqs. ( ^.23| ), we find: 



7? = 37(2- 37)Ao(r!o + Ao-l), 



(5.24) 



where we have used eq. (2.16) and set ao = 1 as usual. 
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Eq. (5.22) is linear as desired, but not straightforward because of its variable coefRcients. We can recast it in normal 

J e~'^^°'^da, where 4>{a) = J{(a^^)da = Clna. This procedure 



form by changing independent variables from a to z ■ 
leads to the following differential equation for y{z): 



where z = (1 — ^ ^ and: 



pi : 



2(2-0 

1-C 



(5.25) 



(5.26) 



Somewhat surprisingly^ eq. ( 5.25| ) has the same form as the differential equation (3.7) governing the solutions of the 
A (X models in §|l|. (The two are identical if we put z ^ t, y x, ^ ^ £, and v — )■ —a.) 

So we could in pr inciple brin g over all the results of §[11, for the cases /i = 1,2,3,4 at least. However, combining 
the definitions ( 5.23| ) and (5.26), we find that = (2/37) (37 - 1) ; or 7 = [3(1 — /i/2)] ^. Therefore solutions obtained 
in this way would correspond to equations of state with 7-values of 2/3, 00, —2/3 and —1/3 respectively. These do not 
describe realistic forms of matter, at least not in the present universe ^T\ . Conversely, values of 7 that are reasonable 
(suc h as 7 = 1 or 4/3) correspond to non-integral values of /x (such as 4/3 and 3/2 respectively). It is doubtful that 
eq. (5.25) can be solved analytically in these cases. We therefore leave the possibility that n = 4 for future numerical 
analysis. 



VI. A AS A FUNCTION OF THE DECELERATION PARAMETER 



A. Evolution of the Scale Factor 



We turn finally to the last of our phenomenological decay laws, eq. (2.13), writing it in the form 



A^V [ - 

a 



(6.1) 



As far as we are aware, no such dependence has previously been considered for the cosmological term. However, it is 
a natural extension of the other decay scenarios considered so far. There is no fundamental difference between the 
first and second derivatives of the scale factor that would preclude the latter from acting as an independent variable 
if the former is acceptable. 



Substituting the decay law (S.l) into eq. (2.£), we find: 



da 



7 \ da 



(6.2) 



We adopt the value r 1 for the remainder of § VI , si nce w e wou ld li ke to solve for the Hubble parameter in analytic 
form in order to make use of of the Taylor expansion (O). Eq. ( |6.2| ) then takes the form: 



dH / 3-V \ k 
da a\jV-2j a^\-fV-2 



37-2,^_,^ 



(6.3) 



As in § V_B, let us make a change of dependent variables from H to x = H'^ , where s is an arbitrary constant. Eq. (6.3) 
then takes the form: 



dx 
da 



57/3-1? 



sk / 37-2 
^ V7X>- 2 



a \jV-2 

If we choose s = 2, this is reduced to the linear form ( |5.6| ), with: 

X>- 3 



X 



{s-2}/s 



S{a) = 27 



7X>- 2 



T{a) = 2k 



37-2 
7X>- 2 



(6.4) 



(6.5) 



22 



Multiplying through by a factor of exp[J S{a)da] 
Hubble parameter: 



Hia) 



(227(15 3)/(fV 2) g^^^ solving exactly as in §VC, we obtain for the 

1/2 



a 



(6.6) 



where Cq is a constant of integration. Imposing the boundary condition H{ao) — H{1) — Hq, we find with the help 
of eq. (2.16) that Co — Hq + k ~ i/Q(f7o + Aq). We can also fix V in terms of observable quantities. With r = 1, 
the decay law (OA) gives D = Ao(a/a)t=tp = — Aq/H ^ qq — — SAp /gp, where qo is the present value of the deceleration 
para mete r, and we have used the definitions ( ^.1E| ) and (2.18). Substituting this result into eq. (6.6) along with 
eq. ( ^.16 ), and recalling that H = {Ho/a)da/dT, we find: 



where: 



da 
d7 



(r!p + Ao) , 



1/2 



a 



(2 - 37)gp 
go + (37/2)Ap ' 



(6.7) 



(6.8) 



This expression, together with its time derivative: 



(6.9) 



can be substituted into the Taylor expansion (4.9). 



B. Minimum Values of the Scale Factor 



As in 
implies: 



VD, we require that oscillating models satisfy da/dr = at a = a,. In conjunction with eq. (6.7), this 

1^0 + Ap - 1~ 



ilp + Ap 



(6.10) 



Let us write this out explicitly using the last of eqs. (|6 
two possible cases: 



As with eq. (5.15), we will find it convenient to distinguish 



^0 + Ap 



S2o + Ao - 1, 

l^p + Aq - 
5ip + Ap 



+ (3-fAo/29o) 1 
3t-2 J 



(37Ao/2go)-1 1 
37-2 J 



if qo > qc 
if qo < qc, 



(6.11) 



where we have defined qc = — (37/2)Ap (leaving the case qo = qc aside for t he time being). 

We can draw a number of useful conclusions from the form of eq. (6.11). Firstly, (1) that spatially flat solutions 
(1 7p + Ap = 1) again have either a,, — oo (if qo > qc) or a* = (if qo < qc)- This is just as in the previous section 

{&■ 

2) Secondly, requiring real values for a* (subject to the same proviso about even-numbered integer exponents as in 



j), we can conclude that closed solutions (ie., Ap > 1 — Qo) must satisfy Ap > — ilp, while open ones (Ap < 1 — flo) 
obey Ap < — fip. It follows that, if Ap is a positive quantity, as observations almost certainly indicate (§IVD2), then 
models with a nonzero minimum scale factor must, in the present theory, be closed. We will assume that both these 



conditions hold in the remainder of S VI 



(3) Thirdly, requiring that < a* < 1 as before, we learn that the deceleration parameter qo satisfies: 

qa < qc- 



(6.12) 



This follows from the fact that jfip + Ap — 1| cannot be greater than |51p + Ap| (assuming that flo and Ap are both 
positive) . 

(4) Finally, combining the conclusion (3) with the definition of qc, we infer that: 
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go < -(37/2)Ao. 



(6.13) 



Assuming as we are that the cosmological term is positive, this imphes that the deceleration parameter must be 
negative for a universe filled with normal matter (7 > 2/3). 

Unfortunately, the deceleration parameter go remains among the most poorly-constrained quantities in observational 
cosmology. Nevertheless, it is fair to say that the majority opinion among cosmologists holds that go is probably 
positive |99| ]. The most recent experimental determination, obtained from Type la supernovae, leads to a value of 
go = 0.385 ± 0.36 [100|. Since oscillating models in the present theory not only have go < 0, but fc > as well, they 
are somewhat disfavored in comparison to those of §IV and §0; and we judge that this is a reasonable place to halt 
our investigation for the time being. 



VII. CONCLUSIONS 



We have examined the evolution of the scale factor a{t) in the presence of a variable cosmological term A, and also 
extended existing treatments by adopting a fairly general equation of state for ordinary matter. 

A number of new exact solutions for a{t) have been obtained in cases where A cx t~^. These models are singular, 
but can be significantly older than those in which A = constant. For odd values of £, the cosmological term must be 
negative (or zero) if the scale factor is to be real-valued. Our conclusions may not extend to cases in which fc 7^ 0; 
this will require more detailed numerical analysis. 

For a cosmological term that scales as A c>c a~™, we have solved numerically for the scale factor as a function of 
time, and found that there are closed models which are compatible with observation and contain no big bang. This is 
in sharp contrast to the situation where A = constant, for which experimental evidence firmly establishes the existence 
of an initial singularity. (The variation effectively allows one to obtain a large A-term where it is most important 
— near the "big bounce" — without the price of a large cosmological constant at present times.) This appears 
not to have been widely appreciated, probably because variable cosmological terms have so far been studied almost 
exclusively in the context of the cosmological "constant" problem. We have obtained constraints from experimental 
upper limits on Ao, as well as requirements of sufficient age, normal gravitational lensing at high redshifts, and others. 
As specific numerical examples, oscillating models with zero pressure, fto = {0.34, 0.51, 0.61, 0.67} and Aq in the ranges 
{0.68 - 0.72, 0.51 - 0.57, 0.41 - 0.49, 0.39 - 0.48} are observationally viable if m = 1, 1.5, 1.8 or 1.9 respectively. If the 
bounce is to be deep enough to generate the temperatures required by conventional nucleosynthesis, then m > 1.8. 

We have also solved numerically for the scale factor when A oc H". In this case we have found open models which 
can account for the observational data despite their lack of an initial singularity. In particular, oscillating models 
with zero pressure and values of Ao either just above or just below 1 — ilo ^-i'^ viable if CIq < 2/3 and n = 2. 

(We have investigated only the cases n = 2 and 4 in detail.) If fig « 0.3 and Aq « 0.15 (or 0.7), for example, then 
the most recent "big bounce" could have been deep enough to account for phenomena such as the cosmic microwave 
background radiation and light element synthesis in a model with n = 2. 

For a cosmological term that depends on the deceleration parameter via A oc g*", we have solved only the case r — 1. 
Closed oscillating models are possible, but require that g be negative at the present time, if the cosmological term is 
positive. 
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FIG. 1. Evolution of the scale factor for flat models with 
A oc r^^ and 7 = 1. Values of Aq are labelled beside each 
curve, and Slo = 1 — Ao in each case. 

FIG. 2. Evolution of the scale factor for flat models with 
A oc r~ and 7 = 1. Values of Ao are labelled beside each 
curve, and Slo = 1 — Ao in each case. 

FIG. 3. Evolution of the scale factor for flat models with 
A oc r"'' and 7 = 1. Values of Ao are labelled beside each 
curve, and Slo = 1 — Ao in each case. 

FIG. 4. Evolution of the scale factor for flat models with 
A oc r"'* and 7 = 1. Values of Ao are labelled beside each 
curve, and fio = 1 — Ao in each case. 

FIG. 5. Evolution of the scale factor for models with 
m = 0,7 = 1, f2o = 0.34, and values of Ac labelled beside 
each curve (after Felten and Isaacman [5]). 

FIG. 6. Phase space diagram showing constraints on mod- 
els with m = and 7 = 1 (after Lahav et al [15]). 

FIG. 7. The age constraint tq > 0.5: (a) as a function of m, 
assuming 7 = 1; and (b) as a function of 7, assuming m = 0. 

FIG. 8. The lensing constraint za > 4.92: (a) as a function 
of m, assuming 7 = 1; and (b) as a function of 7, assuming 
m = 0. 

FIG. 9. Enlarged view of the phase space diagram. Fig. ^ 
now plotted for various values of m between and 2 (labelled 
beside each pair of curves), assuming 7 = 1. 

FIG. 10. Evolution of the scale factor for universes with 
m = 1,7 = 1 and Qo = 0.34. Compare Fig. ^. 

FIG. 11. Phase space diagram for the case n = 2 with 
7 = 1, showing contours of equal minimum size a*. 

FIG. 12. Evolution of the scale factor for models with 
n = 2,7 = 1, f2o = 0.3, and values of Ao labelled beside 
each curve. 
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TABLE I. Examples of phenomenological A-decay laws. 



dX/dt oc A" 



Decay Law^ Reference 

AcxT'' l"] 

AocT" [21] 

Aoce""" [42] 
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A (X a-^(i+'=) 47 
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A (X a"™ 
dA/df oc aiy"A 
dA/df oc //^ 

A cx C + /3a-™ [h8p 



[«1 
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~1 
~1 



Aoct*-2 + /3t2(<-i) [6C] 

A oc /3a-2 + -ff^ [|ll62l 

A oc + ^i-2/« |gj] 

AocC + c-'" [El~l 



A oc C + /3a-2 + _9-2 |6E 

A oc /Sa-"" + \~t] 

A oc IL" 

Aoc{l + l3H){H'' + k/a') m 

Aoct-\i3 + ty in] 

dA/dt oc /3A - A^ "] 

Aoca"^ |~1 

Aoca"^ + /3a-* |~1 

A oc + l3aH{dH/da) |7 ' 



^r, a, are the temperature, scale factor, time and Hubble 
parameter respectively; while /?, e, £, m and C are constants. 



TABLE IL Borner-Ehlers-type upper limits on matter den- 
sity Qo for various values of m and 7, assuming 2, > 6. 

7: 2/3 5/6 1 7/6 4/3 5/3 2 



m: 






00 


0.033 


0.006 


0.002 


0.000 


0.000 


0.000 


1/4 


00 


0.16 


0.10 


0.079 


0.067 


0.053 


0.044 


1/2 


00 


0.31 


0.20 


0.16 


0.14 


0.11 


0.090 


3/4 


00 


0.49 


0.33 


0.26 


0.22 


0.17 


0.14 


1 


00 


0.75 


0.48 


0.38 


0.32 


0.25 


0.20 


5/4 


00 


1.1 


0.71 


0.55 


0.46 


0.35 


0.29 


3/2 


00 


1.9 


1.1 


0.87 


0.71 


0.54 


0.44 


7/4 


00 


4.2 


2.4 


1.8 


1.4 


1.0 


0.84 


2 


00 


00 


00 


00 


00 


00 


00 
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